Congestion phenomena on complex networks 
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We define a minimal model of traffic fiows in complex networks containing the most relevant 
features of real routing schemes, i.e. a trade-off strategy between topological-based and traffic-based 
routing. The resulting collective behavior, obtained analytically for the ensemble of uncorrelated 
networks, is physically very rich and reproduces results recently observed in traffic simulations on 
scale-free networks. We find that traffic control is useless in homogeneous graphs but may improves 
global performance in inhomogeneous networks, enlarging the free-flow region in parameter space. 
Traffic control also introduces non-linear effects and, beyond a critical strength, may trigger the 
appearance of a congested phase in a discontinuous manner. 

PACS numbers: 02.50.Ey, 68.35.Rh, 89.20.Ff, 89.75.Fb 



The first identified Internet's congestion collapse dates 
back to October 1986, when data troughput from 
Lawrence Berkeley National Laboratory to the Univer- 
sity of California in Berkeley dropped from 32 Kbps to 
40 bps. After that initial event, traffic congestion con- 
tinued to threaten Internet's practitioners, even after the 
implementation of congestion control algorithms able to 
recover the system in case of traffic overloads yjj. Com- 
puter scientists have also elaborated several schemes of 
congestion avoidance^ that should prevent congestion to 
occur by keeping the system far from high levels of traffic 
0. Congestion avoidance and control are performed by 
continuously updating the dynamics of end-to-end flows 
in response to the variation of the load level in the net- 
work. Their functioning depend on the average round- 
trip-time (RTT) of the Acknowledgement signals (ACKs) 
used to exchange information between routers. For this 
reason, the observation of heterogeneous patterns in RTT 
time series has been often proposed as evidence of peri- 
ods of congestion 0]. Apart from these indirect measures, 
congestion events are difficult to monitor and study, so 
that a clear phenomenological picture is still missing. In 
spite of the wide interest in developing optimal routing 
algorithms, much less attention is devoted to explore the- 
oretically the dynamical mechanisms responsible of con- 
gestion. It is possible that a better comprehension of 
these mechanisms could help in understanding experi- 
mental data, in predicting congestion events and design- 
ing better routing protocols. 

The topological and dynamical properties of dis- 
tributed information systems, such as the Internet 
pose theoretical challenges of a similar nature of those 
addressed in statistical physics. Therefore, understand- 
ing network congestion phenomena has become a subject 
of intense research in this field Q , in particular after the 
works by Takayasu and collaborators Q, in which the 
evidence of a phase transition from a free-flow regime 
to a congested phase depending on the load level was 
reported. In a recent work, Echcniquc et al. [3] have 



shown using numerical simulations that the nature of 
the congestion transition depends on the type of routing 
rules. They have adopted a routing scheme that could 
be considered a first approximation of realistic transmis- 
sion control protocol (TCP) routing: packets follow the 
shortest path between their source and destination, but 
small detours are admitted in order to avoid congested 
nodes. In case of purely topological routing (e.g. along 
the shortest paths) they found that the congested phase 
appears continuously, whereas the transition is discon- 
tinuous if some traffic-aware scheme is considered. The 
effect of routing rules on network performance has also 
been addressed in 

In this Letter, we put forward a minimal model of traf- 
fic that preserves all interesting features previously ob- 
served in simulations but is simple enough to be studied 
analytically. Both continuous and discontinuous phase 
transitions observed in 7] are reproduced, and their re- 
lation to microscopic packets dynamics is clarified. 

Let us consider a network of nodes and let v{i) de- 
note the set of neighbors of node i. We describe traffic dy- 
namics as a continuous time stochastic process, in which 
packets are generated at each node i with a rate pi. Each 
node is endowed with a first-in first-out (FIFO) queue in 
which packets are stored waiting to be processed. Let 
Ui be the number of packets in the queue of node i. If 
n.i > 0, node i attempts to transmit packets at a rate 
r^, which represents bandwidth, to one of the neighbors 
j € v{i). We assume the following probabilistic routing 
protocol. First, the node j is chosen at random among 
the neighbors v{i) of i. Second, the fate of the packet 
being transmitted depends i) on whether j is the desti- 
nation node for that packet and ii) on the state of conges- 
tion of node j. We model both as probabilistic events: i) 
we call /ij the probability that node j is the destination 
of the packet being sent, meaning that with probabil- 
ity jij the packet is "absorbed" in the transfer; and ii) 
we assume that the transfer is refused by node j with 
a probability rjirij)^ which is increasing with Uj] in this 
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case the packet does not leave node i. This model relates 
in a stylized manner 1) the structural features, encoded 
in the network and in nodes' characteristics (j>i, /i^ and 
Ti), 2) the protocol response to traffic, detailed in the 
function ri{n) and 3) the ensuing traffic process. 

Our model is a simplifyied version of more elaborate 
models0, 0]- In these, packets are created at each node 
with a given destination node, randomly chosen. Then 
they are dispatched by a given routing protocol, minimiz- 
ing some cost along the path. It is reasonable to assume 
that the statistical nature of the collective behavior does 
not depend crucially on the details of the routing proto- 
col, and that it can be captured by the much simpler ran- 
dom diffusion-annihilation process assumed by our mode. 
In particular, the probabilistic traffic-aware routing as- 
sumed above is a simple analytically tractable way to in- 
troduce a traffic dependence in packets' dynamics. Also, 
in shortest-path routing a node of degree ki = \v{i)\ is 
visited with a probability oc fcf , with /3 w 2, whereas 
in the random walk protocol assumed by our model the 
probability of visiting node i is proportional to the de- 
gree ki. In both, high degree nodes are more exposed to 
events of congestion, which is why star-like topological 
structures are particularly vulnerable to congestion and 
perform optimally only at low traffic [§]. Our model, can 
easily accommodate for this statistical features with a 
specific degree dependence of or by considering degree- 
biased random walks, like in [lo| . 

The phase transition from free-flow to a jammed phase, 
is characterized [7.] by the order parameter 

i.e. the percentage of not adsorbed packets for unit time, 
where M{t) = ni{t) is the total number of packets in 
the system at time t, P — J^i Pi is the rate of creation of 
packets and r is the observation time. Note that a local 
order parameter, replacing Af{t) by and P hy pi, 

can be defined in the same way. 

The model discussed above can efficiently be analyzed, 
in the asymptotic time limit, within a mean field approx- 
imation Vijii, . . . , njv) — Y\i T^iini)- The corresponding 
master equations can be solved by message passing-type 
algorithm [ill] . 

Rather than investigating specific examples, we prefer 
here to focus on the mechanisms inducing the change 
from a continuous to a discontinuous phase transition 
when traffic-aware routing is considered. A major insight 
is obtained rephrasing the problem in terms of ensembles 
of graphs. We consider uncorrelated random graphs with 
degree distribution P{k), so that Uk represents now the 
average queue length of nodes in classes of degree k. We 
focus on the simple case iJ,i — fj,, pi — p and Vi = 1 for 
all i, and the routing protocol r]{n) — for n < n* and 
•q{n) = f] for n > n* . The mean-field transition rates for 
nodes with degree k are Wk{n n-\-l) = p + {1 — fi)(l — 




FIG. 1: p{p/fi) for an homogeneous graph from theoretical 
predictions for rj = 0.25, 0.75. Inset: phase diagram for the 
same graph. 

q)!^[l-fj6{n-n*)] and -^71-1) = 6{n){l-x)- Here 
z is the average degree, qk is the probability that a node of 
degree k has empty queue and Xk = vP{ni > n*\ki = k} 
is the probability that a node of degree k refuses packets. 
Likewise q = J2k <lkP{k) and X = J2k ^XkP{k). 
The average queue length Uk follows the rate equation 

hk = p + (1 - /i)(l - g)-(l - Xk) - (1 - «fe)(l - X) (2) 

z 

Note that summing over k and dividing by p we obtain a 
measure of the order parameter p{p). Since hk depends 
linearly on fc, high degree nodes are more likely to be con- 
gested. Therefore, in the stationary state for a given p, 
there exists a threshold real value k* such that all nodes 
with k > k* are congested whereas nodes with degree 
less than k* are not congested. Congested nodes {k > k*) 
have qk — and Xk = V- For k < k* , the generating func- 
tion Gk{s) = J2n P(.''^k)s^ of the packets distribution can 
be computed from detailed balance. This takes the form 

Gkis) = Qk { + corresponding to a 

double exponential, where Ofc ~ x] 
and bk — ^[(1 — ~ ^ x]- From the normal- 

ization G'fc(l) = 1 and the condition rik = 0, we get 
expressions for qf^, Xk and, finally, for q, x- 

The value k* is self-consistently determined imposing 
that nodes with k = k* have g^* = 0, Xk' — V and 
•hk* = that translates into the equation 

fc* = [l-p-x]/[(l-M)(l-^)(l-9)] (3) 

The above set of closed equations can be solved nu- 
merically for any degree distribution P{k) and p{p) can 
be accordingly computed. The solution is particularly 
simple for regular graphs: ki = z, Vi in the limit 
n* ^ 1. Then the congestion free solution with p — 
has qk = q = 1 — and Xfc = X = and it exists 
for p < p. In the congested phase, instead, all nodes 
have n,; — !■ cx), i.e. x = V and q = 0. This solution has 
p = h/p = 1 — (1 — fj)p/p and exists for p > (1 — fj)p. 
Therefore, in the interval p € [(1 — fj)p,p] both a con- 
gested and a free phase coexist, as shown in the inset of 
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uncongested nodes exist up to kp = k*{l — rf). Using 
this classification, we get a first expression for x, i.e. 



Xi = E 
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FIG. 2: p{p) for an uncorrelated scale-free graph {P{k) oc 
A:-^ k^in = 3, kmax = 110, = 3000), fi = 0.2, n* = 10 and 
fj = 0.05 (below, green) and f; — 0.95 (above, red), from both 
simulations and theoretical predictions. Inset: Hysteresis ci- 
cle for the same graph for r; — 0.95. 



Fig. [T] The behavior of p as a function of p exhibits 
hysteresis: The system turns from a free to a congested 
phase discontinuously as p increases a.t p — ^ and it re- 
verts back to the free phase from a congested phase at 
p — {1 — f])^ as p decreases. This simple case also shows 
that traffic control is useless in homogeneous graphs, as 
it does not enlarge the stability region of the free phase, 
while making a congested phase stable for p S [{1—f])^, p] 
(see inset of Fig. [T|). 

The case of heterogeneous graphs instead is much 
richer. In Fig. [2l we display p{p) for a scale-free net- 
work from both simulations (points) and numerical cal- 
culations (full line) . The agreement is very good and the 
behavior of the curves reproduces the scenario observed 
in The figure is obtained for = 0.2 and n* — 10, 
but the behavior does not qualitatively change for dif- 
ferent values of these parameters. The dependence on 
fj brings instead qualitative changes. Increasing f] from 
0.05 to 0.95, the transition becomes discontinuous and 
Pc increases. FigI5] (inset) also shows that in case of 
discontinuous transition, the system exhibits hysteresis 
phenomena: a congested system does not immediately 
decongest if the creation rate p is decreased under the 
threshold value Pc- 

When the whole system is congested, simple arguments 
of queueing theory show that p{p) follows the curve 
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; however the most interesting situation occurs 
when the network is only partially congested. This case 
can be better understood considering the limit n* — > oo, 
in which the calculations simplify considerably without 
modifying the overall qualitative behavior. 
In this limit, uncongested nodes have Uk < 1, hence 
Xk ^ and qk = 1 — a^, meaning that these nodes 
have short queues and do not reject arriving packets. If 
afe > 1, then Qk ^ and Xk — ^ "^"^ ; more precisely 

we have Xk = ^ ^ iri^ ^ ^) k < k* and Xk ^ fj 
for k > k* . The latter class identifies congested nodes, 
while we call fickle those with k*{l — fj) < k < k* . The 
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iPik) +fiEtr^ iPik). Eq. 
^ provides a further relation between q, x and k* . q 
can be eliminated using its definition wich leaves us with 
an implicit equation for whose solution we call X2- 
In Fig. Owe plot the difference Ax = Xi — X2 vs. k* , 
for fj = 0.05 (left) and 0.95 (right) and different values 
of p. The zeros Ax(fc*) correspond to the only possible 
values assumed by fc*. In both cases, Ax decreases as 
p. For small rejection probability (77 = 0.05 in Fig. [3]), 
there is only one solution k*(p), which decreases from 
+00 when increasing p from 0. The value Pc at which 
k*{pc) = kmax is the critical creation rate at which high- 
est degree nodes become congested. At larger p, k*{p) de- 
creases monotonously until all nodes are congested when 
k*{p) = kmin- Hence for low values of fj, the transition 
from free-flow to congested phase occurs continuously at 
the value of p for which k* (p) — kmax ■ 
At large fj {fj — 0.95 in Fig. [3]), the scenario is much 
more complex. Depending on p, the equation can have 
up to three solutions, kl{p) < k2{p) < k^{p). It is easy 
to check that only fcj and fcg are stable solutions. For 
p <C 1 there is only one solution at k^{p) ^ k^ax, cor- 
responding to the free phase. As p increases, two cases 
are possible: i) k^{p) reaches kmax before the local min- 
imum crosses zero, in which case the congested phase 
emerges continuously; ii) if instead the solution kKp) 
appears when still k^(p) > kmax, a congested phase ap- 
pears abruptly. For sufficiently large values of p, only 
kl{p) survives, and the network becomes fully congested 
when kl{p) < kmin- Therefore, the existence of a purely 
discontinuous transition depends strongly on the tail of 
the degree distribution. 

In the latter case ii) the hysteresis phenomenon becomes 
evident upon varying p in opposite directions: Starting 
from the free phase at low p, the system selects the solu- 
tion fcg (p) and follows it upon increasing p until the so- 
lution k^{p) disappears. On the contrary, starting from 
the congested phase (large p) the system locks into the 
congested phase kl(p) and remains congested until the 
solution kl{p) disappears, with a discontinuous transi- 
tion (see inset of Fig. [2]). Hysteresis also occurs in case 
i) for the same reasons. A detailed account of this rich 
phenomenology will be reported elsewhere jlli] . 

We computed the phase diagram (see Fig. |4|) in the 
plane {p, fj) for the same uncorrelated scale-free random 
networks with 7 = 8 considered in Fig. [21 in the case 
n* —^ 00. The dashed line represents the continuous 
phase transition, separating free-fiow regime from con- 
gestion. At the point C, the critical line splits in two 
branches that define a coexistence region. The upper full 
line represents the discontinuous transition from the free- 
phase to the jammed state, whereas the lower indicates 
the opposite transition from congestion back to free-flow. 
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FIG. 3: The zeros of Ax(p) vs. k* define the threshold degree 
for the onset of congestion in a network. The picture refers 
to a scale-free random network with 7 — 3.0 {kmax = HO), 
and different values fj = 0.05,0.95 and p. The solution kl{p) 
in the right panel falls outside the plot. 
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FIG. 4: (?7,p) phase diagram for an uncorrelated scale-free 
graph (P(fc) oc fc-^ kmir. = 3, k^ax = 110, N = 3000), 
fj, = 0.2, n* = 00 from theoretical predictions. The inset 
shows PciN) for 77 = 1 (red line) and rj = (black line). 



The dotted line decreasing from the maximum of the crit- 
ical Une, is an unphysical branch of the solution obtained 
from calculations. Indeed if a free flow phase is stable 
for a given value of fj, this will persist for larger values 
of fj because traffic control only affects congested nodes 
in the limit n* 00 (i.e. a free stationary state cannot 
become congested if we increase fj) . 

This phenomenology crucially depends on the tail of 
the degree distribution. Since kmax depends on the sys- 
tem's size, we expect that Pc depends on N as well; the 
inset of Figd] show that for the same scale-free network 
of Fig. [3] the critical rate of packets' creation goes to zero 
as Pc{N) cx 1/^/N, for 77 = 0, but it goes to a constant 
for ?7 = 1, in the limit of large N. Hence the point C 
separates two regions in fj with distinct behavior of finite 
size effects. 

The mechanisms triggering the emergence of conges- 
tion is somewhat reminiscent of jamming or bootstrap 
percolation, where a node is occupied if the number of 
occupied neighbors exceeds a given threshold. In these 
models, as the threshold increases, the transition turns 
from continuous to discontinuous ■ Similarly, as traf- 



fic control in the routing protocol is enhanced, the con- 
gestion transition may turn from continuous to discon- 
tinuous, as a result of a cooperative process in which the 
more congested neighbors a node has, the more it is likely 
that it will be congested. 

In conclusion, we have proposed a minimal model to 
study the emergence of congestion in information net- 
works. For un-correlated random graphs, the analy- 
sis can be performed analytically at the ensemble level, 
revealing that the interplay between the feedback pro- 
cess induced by traffic-aware routing and the topologi- 
cal structure of the network (in the tail of the degree 
distribution)_generates a rich phenomenology of phase 



transitions 
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Traffic-aware routing is useful only in 
heterogeneous networks, where it expands the region of 
stability of the congestion free state. However, when its 
effects are strong enough a congested phase may arise 
abruptly, and once it arises, it may persists even under 
lower traffic loads. 
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